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Abstract 

Form factors and decay widths for D* Dy and D* —> Dn decays are 
estimated in a relativistic constituent quark model. Relativistic corrections 
due to light quarks are found to be substantial and to suppress the vector and 
axial vector form factors. The CLEO experimental value of = r(D*^ —> 

= 0.572 ± 0.057 ± 0.081 is used to determine the quark 
axial vector coupling qa-, which is found to be 0.6—0.8 for = (350—200)MeU 
correspondingly, as compared with the chiral model result gA = 0.8 — 0.9. The 
heavy meson-pion strong coupling g is found to be 0.4 — 0.6, much smaller than 
5 = 1 which is expected in the large Nc and nonrelativistic limit, but consistent 
with some heavy hadron chiral theory and QCD sum rule results. 


Because the chiral SU{2) symmetry is spontaneously breaking, the axial vector 
coupling of the constituent quark, qa, may take any values. It has recently been 
argued based on a chiral model^^l that to leading order in 1/Nc (where Nc is the 
number of colors in quantum chromodynamics ), the constituent quarks just behave 
like bare Dirac particles, and both u and d quarks have vanishing anomalous magnetic 
moment and have the axial vector coupling 

9A = 1. (1) 

With order I/Nq corrections, the value of qa is slightly reduced and is estimated to 
be[2-3] 

gA = 0.8 - 0.9. (2) 

On the other hand, using the experimental value for the nucleon axial vector coupling 
Ga = 1-25 with the static SU(6) wave function for the nucleon leads to qa = 1.25x | = 
0.75. The discrepancy between this value and (2) might be removed by including the 
relativistic corrections to the nucleon wave functions^. 

The heavy meson decays, e.g., D* ^ D decays may also be a good testing ground 
for the value of quark axial vector coupling qa- Although the total widths of D*~^ and 
D*^ are still unknown, their decay branching ratios have been measured accurately 
by the CLEO Collaborationt^l (see also ref. [6]). These heavy meson decays have been 
discussed in the hterature['^“33]_ 'Jq determine qa, the most useful data might be 
B{D*^ — D^'j) and B{D*^ which are better measured experimentally and 

are not sensitive to the value of the charm quark mass in theoretical calculations. We 
suggest using the measured ratio^^] 

= 0.057 ±0,081, (3) 

together with the calculated r{D*^ D^l) and the form factors in the hadronic 
decay to determine qa- 

Because the light constituent quarks {u and d) inside the heavy mesons are rela¬ 
tivistic, to discuss the qa issue in the heavy meson decays the relativistic motion of 
light quarks has to be taken into consideration. In the following we will study the 
heavy meson decays and the qa problem in a relativistic constituent quark model, 
based on the Bethe-Salpeter (BS) formalism. 

For a qiq2 bound system with quark momentum gi, mass mi; antiquark momentum 
q 2 , mass m 2 ; relative momentum q; and total momentum P, meson mass M; and 

7Tl • 

Pi = rjiP + q, p2 = P2P -q, Pi = - y - (* = 1,2), 

mi -|- m2 

the general form of the three dimensional BS wave functions for the 0“ and 1“ mesons 
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is given by (see e.g. ref. [14]) 


4>p(^ = A^pi)7°(1 +^)757°A^(P2)0(^ (4) 

4>p(^ = A^p 1)7°(1 +^)e7°A2(p2)/(^ (5) 

where is the polarization vector of the 1 “ meson, = 0 ; i?* = \jpi^ + mf, and 

A!|_ = |(1 ± (iP+mi) ^ positive (+) and negative (-) energy projectors, and 0, 

/ are scalar wave functions. (4) and (5) respect the space reflection symmetry, and 
also respect the flavor-spin symmetry in the heavy quark limit. This can be seen by 
taking mi oo, then pi and (4) and (5) become 

+^)^57°A^(P2)0(^), ( 6 ) 

v'j 

= ^(1+^)e7°A^(pl)/('?), (7) 

where = ^, and f = <p which is due to the vanishing of the color magnetic force 
in the heavy quark limit. The normalization condition for (4) and (5) is 

^ f d^qTr = 2P° = 2\/P^ + M^. (8) 

zn J L 

Suppose a flavor changing quark operator g/Tg* induces the transition 

^p{^{pi,ml-,p2,m2]P,M) ^ ^p>{q'){p[,m[-,p2,m2]P',M'), ( 9 ) 

where the antiquark {p 2 , m 2 ) remains a spectator, then the transition matrix element 
is given by 

(P' I q-fVq, \P) = ^J Pp^Tr [4>l„(g7)7°r4>p(^] . (10) 

Likewise, the matrix element induced by the antiquark transition can be written in a 
similar way. 

In the radiative Ml decay such as D* —> P 7 if we neglect the quark anomolous 
magnetic momentl^’^l we then only need to calculate the matrix element induced by 
the vector operator gy^g of the quark and antiquark, and according to (10), (4), and 
(5), for the quark transition we have 


= {P' I I P^ e) 


/■ ,3 rj. P2 -m2 . p.p'i + mi Pi +mi 

2 E 2 75 ( 1 +^) 2E[ 2Ei 


^ m2 B\z / ^ "^2 Bo 

xfi{P 2 -^- P)<Pf{P 2 -^- P) 

mi + m2 mi -|- m2 
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where v = -^iV' = In general, the vector current matrix element (11) can be 
expressed as 

jf, = , ( 12 ) 

where is the vector form factor due to the quark qi transition. If the scalar wave 
functions fi (for the initial state 1“ meson) and 0/ (for the dual state 0“ meson) in 
(11), which are to be determined by the interquark forces, are known, then a direct 
calculation for (11) will give the form factor ,^yi. We will give this result later on by 
solving the BS equation with a QCD-motivated interquark potential. Before doing 
that calculation we may first consider a simplified calculation as follows. Because in 
the heavy quark limit the color-magnetic force vanishes, the 0“ and 1“ mesons will 
have the same spatial wave functions, we may assume fi and to take the same 
form 

fi{p) = aif{p), <Pf{p) = aff{p), (13) 

where a* and a/ are normalization factors determined by (8) in the initial and final 
meson frames respectively. Although in (11) the light quark is relativistic, in order to 
see the relativistic effects more explicitly (but less rigorously) it might be instructive 
to make a nonrelativistic expansion in terms of the inverse of the quark masses. Then 
to the first order we find 
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8mim2 



1 

6 


m2 

mi + m2 
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(14) 


where K = P' — P is the recoil momentum, {(f) = J d^q\ f{f ff is the mean value 
of the quark momentum squared, and (r^) is the mean value of the radius squared 
of the mesons. From (14) we may find some qualitative feature of the relativistic 
effects. We see that the Ml transition can be substantially suppressed if the quark 
is a light quark and becomes relativistic = 0(1). We see also that the effect of 
the spectator antiquark on the quark transition is not strong in general, because the 
coefficient of the term involving {mim 2 )~^ is small. The last term in (14) is due to the 
nonzero recoil momentum, and its contribution is also small because of the smallness 
of the recoil momentum in the D* radiative decays. If the quark is a heavy quark 
then in the heavy quark limit mi oo we will have .^yi = 1 in (14) and 1 and 

therefore (12) will return to the well known expression for the vector current in the 
heavy quark limit. 

In the rest frame of the 1“ meson, F = 0,v° = l,v' = then (12) gives the 
familiar expression for the Ml transition amplitude. With a similar expression for 
the antiquark transition, we can get the 1“ —0“7 decay width for a qQ system 


r, = F? A + 


K 


m. 


rriQ 


(15) 
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where and are the vector form factors for q = u, d and Q = c in the 
D* —> D'j decay. We will calculate these form factors from (11) and will not use their 
nonrelativistic expansion (14), because the light quark can be highly relativistic. 

The hadronic decay D* Dtt can be described by the quark-pion vertex in the 
chiral quark model^^^^ 

(16) 

where qa is the quark axial vector coupling and Tk = 132MeV is the pion decay 
constant. In general, for a qiq 2 system the transition induced by the quark axial 
vector current can be written as 


is 


= {P' I gi 7 ^ 75 gi I P, e) 


(17) 


^Jd^P^Tr 


P2-m2 . -,.Pi + mi Pi + mi 

+ V). 

yfi{p2 - ’^p)Mp2 - ^^P') 

mi + m2 mi + m 2 

iJMd*oMdo {^ai(1 + V ■ v')e^ - ^^2(6 ■ v')v^ - ^A3(e ■ v')v'^} . 


As in the case of .^yi, we can calculate .^ai, 'Ca 2 , and .^a 3 with the scalar wave functions 
/ and (f) obtained by solving the BS equation, but it is useful to give the nonrelativistic 
reduction form for, e.g., the .^ai 


■Cai 


Smf 6 mi + m 2 


(18) 


Again, ,^ai is suppressed by the relativistic motion of the light quark, but the sup¬ 
pression is less severe than the Ml transition form factor ,^yi in (14). This can be 
seen by noting that the coefficient of the term, which is the leading term for the 
suppression, is —1/3 in (18) whereas is —2/3 in (14). 

For e. g. the D*^ D^tt^ decay, | iA |= 44MeI/, | n' |= 0.024, | e ■ v' [-C 1, 
therefore in (17) the contribution of .^a 2 and .^a 3 terms can be neglected, and we then 
get 


r(T)*° ^ DV) 


^A^Al I p ,3 
127rF2 ' " ' 


(19) 


For D*~^ D~^7r^ the width takes the same form as (19), while for D*^ —> D^7i~^ an 

additional factor of 2 on the right hand side is needed due to the isospin difference. 


To calculate the form factors in radiative decay (11) and (15) and in hadronic 
decay (17) and (19), as a simple and naive choice, we hrst use the Gaussian wave 
functions 




( 20 ) 
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for fi and (j)f, where N is the normalization factor, = |(^). For the D and D* 
mesons most estimates give (see, e.g., ref. [16]) 

(^) = 0.2 - 0.3 GeV^ (21) 

Here we will simply take {(f) = 0.21 (0.30) GeV"^ for rUu = rud = 200 (300) MeV, and 
rric = 1500 MeV, and calculate the form factors using their relativistic expressions 
(11) and (17) (not (14) and (18)). The calculated form factors and the 

decay widths are shown in Table 1. We see these numerical results are qualitatively 
consistent with the nonrelativistic reduction expressions (14) and (18) (e.g., is 
more suppressed than .^ai)- Here is determined by using the calculated width for 
D*o ]jo^ g^nd the experimental value for the ratio ~ 0.572 ± 

0.057 ± 0.08115]. 

To be more closely connected with QCD dynamics, we have also calculated these 
form factors and decay widths based on the BS equation with a QCD-motivated 
interquark potentially^!. In the instantaneous approximation the BS equation 

(^1 - mi)xp{q){r^2 + "^2) = ^ y d‘^k G{P, q - k)xp{k), (22) 

where G{P, q — k) represents an “instantaneous” interquark potential in momentum 
space, can be reduced to the following equation for the three dimensional BS wave 
function 

= Jdq\p{q^,q), (23) 

(P° - Pi - = Alf I d^k G{P, q - k)^-^{kh^A^, (24) 

where the contribution of negative energy projectors (i.e. the pair terms) are ne¬ 
glected. 

The interquark potential is described by a long-ranged linear conhning potential ( 
Lorentz scalar Vs ) plus a short-ranged one gluon exchange potential ( Lorentz vector 
By ), i.e. 


V{r)=Vs{r)+^^®rVv{r), 


Vs{r) = \r 


(1 


ar 


Vv(r) = 

3 r 


(25) 


where the introduction of the factor e““” is to regulate the infrared (IR) divergence 
and also to incorporate the color screening effects of the dynamical light quark pairs 
on the QQ potential. It is clear that when ar 1 the potentials given here become 
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identical with the standard linear plus Coulomb potential. In momentum space the 
potentials are 


G(p) = G5(P)+7M®7^Gy(P), 

Gs(p) =--S^(P) + - ^ 

a TT 


{P +a2)2 


3ir"p 


(26) 


where as{P) is the well known running coupling constant and is assumed to become 

^2 

a constant of 0(1) as P —0 


as{P) = 


127r 


27 


ln(a + 


(27) 


QCD 


The constants A, a, a and Aqcd are the parameters that characterize the potential. In 
the computation we will use A = O.lSGeC^, a = O.OdGeV, a = e = 2.7183, Aqcd = 
O.lSOeC. 

Substituting (4), (5), and (26), (27) into the reduced BS equation (24), with quark 
masses = 200 — 350 MeV, rric = 1500 MeV, we can solve for the scalar 

wave functions 0 and / of the 0“ and 1“ mesons respectively, which are different 
due to the color magnetic force induced by one gluon exchange potential. We then 
substitute the obtained (j) and / into (11) and (17) to calculate the form factors. The 
results are also shown in Table 1. 

The form factors are more suppressed for the wave functions by solving BS equa¬ 
tion than for the Gaussian wave functions. This is mainly because, due to the rela¬ 
tivistic correction of one gluon exchange potential (like the Breit-Fermi Hamiltonian) 
in the BS equation the spatial wave function of the 0“ meson has a larger {^) (caused 
by an attractive spin-spin force between the quark and antiquark), and the 1“ meson 
wave function becomes different from the 0“ meson wave function, and therefore the 
overlap integral of wave functions between D* and D is reduced. From Table 1 we 
see that although the decay widths are predicted to be somewhat different in the two 
models, the obtained values for gA_ are close to each other. E.g., for niu = 300MeH, 
gA = 0.65 — 0.67; for ruu = 200MeH, gA = 0.81 — 0.83. Nevertheless, we prefer the 
results obtained from BS equation with QCD-motivated potentials. 

In the heavy quark limit {mq oo) for the hadronic decay D* Dtt the heavy 
meson-pion strong coupling Gd*Dtt, dehned by 


(T>°(p/)7r+(p^) I T>*+(pi,e)) = GD*DnP^ef„ 


(28) 
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may be written in a more convenient form as 


G 


2M, 


D* 


D*Dir — 


-9- 


(29) 


Comparing (28),(29), with (19) it is easy to find 

9 = 9aUi- (30) 

In the large Nc limit ( gA=^ ) and the nonrelativistic limit (,^^i=l ) we would expect 

9 = 1. (31) 


However, based on the relativistic description for heavy meson decays, with a typical 
value of the constituent quark mass = 350(300)MeH our BS model gives 

= 0.67(0.64), gA = 0.60(0.65), g = 0.40(0.42). (32) 

We see that our estimated values of gA are signihcantly smaller than (2): gA = 
0.8 — 0.9, which is expected in the chiral lagrangian approach. This is similar to the 
result of ref. [12]. Moreover, our value for the meson strong coupling g ^ 0.4 is also 
much smaller than (31): g = 1- 

However, our result for the quark axial vector coupling gA is sensitive to the value 
of the constituent quark mass. We see from Table 1 that if = 200MeI/ then the 
BS model calculation would give = 0.58, gA = 0.81. In this case, although the 
light quark inside the heavy meson becomes even more relativistic, the Ml transition 
D*o width gets enhanced due to a larger Dirac moment of the even lighter 

quark. Consequently, with the CLEO ratio (3) the strong decay width of D*^ meson 
and the value of gA will be increased. This possibility is of course not excluded and 
worth further investigating. In this connection it might be interesting to notice that 
rriu ~ 200 MeV is also suggested in some relativistic quark models (see, e.g., ref.[4]). 

As for the meson strong coupling g, with rUu = 200MeH we get g = 0.47, which 
is still much smaller than g = 1. In fact, our BS model calculations show that 
with a wide range of the constituent quark mass, say, = 200 — 350MeV, the 
obtained meson strong coupling is g = 0.47 — 0.40. The Gaussian wave functions give 
g = 0.53 — 0.45 for = 200 — SOOMeV (see Table 1). All these results suggest 

c/ = 0.4-0.6, (33) 


which is much smaller than g = 1 expected in the large Nc and nonrelativistic limit, 
and than g = 0.8 — 1 suggested e.g. in refs.[8,19], but is consistent with the result 
g ^ 0.6 in a heavy hadron chiral models®]. It is interesting to note that some recent 
QCD sum rule analyses favor an even smaller value g = 0.2 — 0.4^7. The study of 
semileptonic decay D ^ ttIui using a chiral effective theory in the heavy quark limit 
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also favors g = Moreover, g = ^ is suggested by some calculation for the 

relativistic effects^^^l. With (29) and (33) we obtain the following estimate for the 
effective D*Dtt coupling and B*Bt[ coupling 


Gd*d-k — 12 — 18, Gb*btt — 32 — 48. 


(34) 


Finally, one might ask whether the nonrelativistic treatment for the light quarks is 
still tenable for heavy meson decays. In the nonrelativistic limit we would have = 
= 'Cai = 1, then from (15) and the CLEO ratio (3) with typical quark masses rric = 
1500 MeV, rriu = rrid = 330 MeV we would get = 107 KeV, = 

IblKeV, which already exceeds the observed upper bound rtoi(F)*’'')(131 
This seems to rule out the possibility that the light constituent quark inside the heavy 
meson can be treated as nonrelativistic, indicating that a relativistic description for 
the light quarks in heavy meson decays is necessary. It is amusing that by above 
nonrelativistic treatment we would get 


gA = 0.72, 


(35) 


which is almost in coincidence with the nonrelativistic value gA = 0.75 obtained from 
the nucleon /5-decay. 

To sum up, we have estimated the form factors and decay rates for the D* —> Dy 
and D* —^ Dn decays, based on a relativistic constituent quark model with QCD 
inspired interquark potentials. Relativistic effects are found to be substantial on the 
form factors. Using the CLEO experimental value 77° = y{I^*o^do^o\ = 0.572±0.057± 
0.081 as input we hnd that the quark axial vector coupling might be consistent with 
gA = 0.8 —0.9, which is expected in the chiral lagrangian approach. However, this can 
only be achieved with a smaller light quark mass, say, rUu = md = (200 — 220)MeV, 
which give larger Ml transition widths accordingly. With the typical value = 
nid = 300 — 350MeU, our obtained value for gA is smaller, say 0.67 — 0.60. It is 
therefore interesting to note that the gA issue might be possibly related to the value 
of the constituent quark mass of light quarks in the heavy meson decays. As for the 
heavy meson-pion strong coupling g, with a wide range of the constituent quark mass 
rriu = fUd = (200 — 350)MeU, we obtain g = 0.6 — 0.4, which is much smaller than 
g = 1 expected in the large Nc and nonrelativistic limit. 




Gaussian 


BS 



ufiu = 200 

niu = 300 

niu = 200 

niu = 300 

niu = 350 


{MeV) 

{MeV) 

{MeV) 

{MeV) 

{MeV) 


0.42 

0.51 

0.36 

0.46 

0.51 


0.97 

0.95 

0.90 

0.87 

0.86 

Cai 

0.64 

0.70 

0.58 

0.64 

0.67 

9a 

0.83 

0.65 

0.81 

0.65 

0.60 

9 

0.53 

0.45 

0.47 

0.42 

0.40 

P(L)*o DO 7 ) 

21 (A:eR) 

15(A:eR) 

16(A:eR) 

13{KeV) 

12{KeV) 

T{D*^ DV) 

37 

26 

29 

22 

21 

Ttot{D*^) 

58 

41 

45 

35 

33 

T{D*+ D+ 7 ) 

0.42 

0.12 

0.24 

0.10 

0.07 

P(D*+ ^ D+ttO) 

26 

18 

20 

15 

14 

t\d*+ D%+) 

55 

40 

43 

33 

31 


81 

58 

63 

48 

46 


Table 1: Form factors and decay widths for D* D'y and D* Dtt decays. Pre¬ 
dicted values are given with (1) the Gaussian wave functions and (2) the wave func¬ 
tions by solving BS equation. Here nic = ISOOMeH, = rud = 200, 300, 350 MeV 
are assumed. The CLEO experimental value of ~ 0-572 ± 0.057 ± 

0.081 is used as input. 
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